Theory of dark resonances for alkali vapors in a buffer-gas cell 



O 



V. Taichenachev,i'2'3 V. I. Yudin,!'^'^ R. Wynands,'*^^ M. Stabler," J. Kitcbing,^ and L. HoUberg^ 

^Novosibirsk State University, Pirogova 2, Novosibirsk 630090, Russia 
^Institute of Laser Physics SD RAS, Novosibirsk 630090, Russia 
'^Time and Frequency Division, NIST Boulder, 325 S. Broadway, Boulder, CO 80305, USA 
^Institut fiir Angewandte Physik, Universitdt Bonn, Wegelerstrafle 8, D-53115 Bonn, Germany 
^Present address: Departement de Physique, Universite de Fribourg, 



Chemin du Musee 3, 1700 Fribourg, Switzerland 
(Dated: February 2, 2008) 



m 
o 
o 

We develop an analytical theory of dark resonances that accounts for the full atomic-level struc- 
' ture, as well as all field-induced effects such as coherence preparation, optical pumping, ac Stark 

, shifts, and power broadening. The analysis uses a model based on relaxation constants that as- 

sumes the total coUisional depolarization of the excited state. A good qualitative agreement with 
' experiments for Cs in Ne is obtained. 

(N ■ 



PACS numbers: 42.50.Gy, 32.70.Jz, 32.80.Bx, 33.70.Jg 



I. INTRODUCTION 



Nonlinear interference effects connected with the atomic ground state coherence are now well known and widely used 
0. One of the most promising classes of these effects, especially for precise measurements, is that of super-narrow dark 
resonances H i3, ,4J that appear in the medium's response to bichromatic laser excitation, when the laser frequency 
^ [ difference is close to the atomic ground-state splitting. The use of vapor cells containing a buffer gas in addition to an 
• ^ ■ alkali vapor has allowed the measurement of resonance linewidths less than 50 Hz 0, While such resonances have 
iCr^' been extensively investigated experimentally (especially in the case of Cs) 0, a detailed theoretical understanding 
(-H , is not yet well developed for realistic multilevel systems, motivating the present work. Our theory was developed in 
O • close connection with ongoing efforts to construct compact atomic clocks 0, Q, H, |1| and magnetometers 0, • For 
any practical application of dark resonances, the stability and accuracy are optimized with respect to parameters such 
as the output signal amplitude, width, and shift. In the problem considered here, many parameters, such as laser 
detunings, field component polarizations and amplitudes, and buffer gas pressure, affect the dark resonance itself. In 
addition, various excitation schemes (for example, D2 versus Di line excitation [lol |) and different atomic isotopes 
' can be used. A natural question arises: what design will optimize the performance of the clock (or magnetometer)? 
, Previous theories did not completely answer this question. One main obstacle was connected with the complicated 

• energy-level structure of the real atomic systems used in experiments. 

I Generally speaking, there are several types of problems in the theoretical description of dark resonances. One 

• problem relates to a proper treatment of the relaxation processes in the system, including velocity-changing collisions 
lll | and the spatial diffusion of coherently prepared atoms 0| . Light propagation through coherently prepared 



^ I nonlinear media, especially through optically thick media |l4j , can be thought of as another type of difficulty. This 
. ^ paper addresses another important problem: that of field-induced processes in multilevel systems such as coherence 
preparation, optical pumping, ac Stark shifts, and power broadening. All existing theories can be classified into three 
kinds: few-state models (basically, three-state lambda systems) |^ Q| i perturbation theories : and numerical 
O ^, simulations All three classes of theories have disadvantages. The first theory neglects many details of the 

• actual configuration of atomic levels. Perturbation theory neglects some effects induced by the presence of the optical 
^ field (namely, optical pumping, ac Stark shifts, and power broadening). Numerical simulation theories demonstrate a 
^ . lack of genuine understanding and predictive power. 

^ ' This paper presents a new analytical theory, that accounts for the level structure (both Zeeman and hyperfine) 
. 5^ , of a real atom, as well as all field-induced effects. The relaxation processes are treated in the simplest way: by 
neglecting velocity-changing collisions and all effects connected with the spatial inhomogeneity, we reduce the model 
to one described simply by relaxation constants. The crucial assumption is total coUisional depolarization of the 
excited state. In addition, we add the (optional) approximations of homogeneous broadening and low saturation. 
With these approximations, a general analytical result is obtained for the atomic response, which result is valid 
for arbitrary excitation schemes {D2 as well as Di lines), light field polarizations, and magnetic fields. In the 
specific case of circularly polarized light in the presence of a magnetic field, where only two states participate in 
the coherence preparation, analytical lineshapes (generalized Lorentzian) coincide exactly with the phenomenological 
model heuristically introduced previously to fit experimental data [l8|. In the case of zero magnetic field, and 
when contributions of different Zeeman sub-states are well overlapped, the resonance lineshape is also approximately 
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is 



FIG. 1: Excitation scheme 

described by the generalized Lorcntzian. A comparison of analytically calculated coefficients of the Lorentz-Lorenz 
model (with no free parameters) with coefhcients extracted from experimental data demonstrates a good qualitative 
agreement. 

II. STATEMENT OF THE PROBLEM 

In this section, the general framework of the problem is described, the basic assumptions we make are stated and 
the specific procedure for calculating the quantities of interest is outlined. We consider the resonant interaction of 
alkali atoms in the 51/2 ground state with a two- frequency laser field 

E(z, t) = El exp[--i(wit - kiz)] + E2 exp[~i{uj2t - k2z)] + c.c. , (1) 

where both components propagate in the positive direction (A:i,2 > 0). The field can excite atoms cither to the P1/2 
state {Di line) or to the P3/2 state {D2 line). Two hyperfine (HF) components are present in the ground state with 
the total angular momenta Fi — 1+1/2 and F2 — 1—1/2 (where I is the nuclear spin). The HF splitting in the 
ground state A = {£1 — £2)/^ is in the range 1 to 10 GHz. The excited state has two {Di line) or four (D2 line) HF 
levels with the angular momenta F^. = I — J^, . . . , I + Je and the energies £e — h-^e- The HF splitting of the excited 
state is typically one order of magnitude smaller than A. To be more specific, wc assume that the frequency lui is close 
to resonance with the Fi — > F^ transitions, while the other frequency uj2 is close to the frequencies of the F2 — > F^ 
transitions. Thus, we have a A-type excitation scheme (Fig.^. In the absence of an external B-field, the HF levels 
are degenerate with respect to the total angular momentum projections. For the Zeeman sub-states the following 
shorthand notations will be used: |e) = \Fe, rrie) with mg — —F^, . . . , F^, and \i, m) = jF^, m) with m = —Fi, . . . , Fi 
(* = 1,2). 

For simplicity, we consider first an atom at rest, positioned at the origin z — 0. Each frequency component of the 
field can in principal induce transitions from both ground-state HF levels. Then the interaction Hamiltonian in the 
dipole approximation contains contributions of two kinds: 

i?Z3-B = - ^ |e)(e|(d.E,)|z,m)(z,m|- ^ |e)(e|(d • EOIj, m)(j, m|e-'("--^)* + /i.e. , (2) 

where we use a rotating frame (the unitary transformation of the ground-state basis \i, m) — > exp(zct;it)|j, to)), and d is 
the dipole moment operator. The first term in ^ is independent of time in the rotating basis, and we refer to it as the 
resonant contribution. The second term, oscillating at the difference frequency, results in off-resonant contributions 
to the optical shifts and optical pumping rates, as well as in temporal oscillations of the atomic density matrix. The 
role of the off-resonant term in the case of a three-level A system has been studied in great detail j3| • The amplitudes 
of the oscillating parts of the density matrix can be approximated as |c?i?p/(?iA)^. For the moderate field intensities 
considered here (< lOmW/cm^) this ratio is very small, \dE\'^ / {fi/S.)'^ ^ 10~® — 10~^, and the oscillating terms can 
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be safely neglected. However, the off-resonant contributions to the optical energy shifts and widths can be significant, 
especially in the case of large one-photon detunings. 

The Hamiltonian for a free atom in the rotating frame can be written as 

Ho = -J2 HSl - ^e)|e)(e| - h^-f m) (1, m| - |2, m)(2, m|) . (3) 

e m 

Here Sj^ = {Si + 52)/2 is the average one-photon detuning, and We are measured from a common zero level (for 
example, from the HF level with maximal momentum — I + Jg), and Sn = S2 — Si = ~ '^i ~ ^ is the Raman 
(two-photon) detuning. 

Since this paper is concerned with the field-induced effects in multi-level atomic systems, the relaxation processes 
are modeled by several constants. The homogeneous broadening of the optical line, due mainly to collisions with a 
buffer gas, is described by the constant 7. We assume that the excited state is completely depolarized due to collisions 
during the radiative lifetime Te, i.e., the depolarization rates 7k obey the condition 

7k T-e > 1 . (4) 

The relaxation of the ground-state density matrix to the isotropic equilibrium, both due to the diffusion through the 
laser beam and due to collisions, is modeled by a single constant T. 

Under the assumption of moderate field intensities and high buffer-gas pressure, we develop the theory in the 
low-saturation limit: 

The two-photon dark resonance appears when the Raman detuning Sji is scanned around zero. The width of the dark 
resonance, which is related to the ground-state relaxation, is usually six orders of magnitude smaller than the optical 
lincwidth 7. The approximation (5ij <C 7 is therefore suitable. 

It should be stressed that all approximations are well justified for typical experimental conditions. For example, 
in the case of Cs in a background Ne atmosphere at a pressure of p — 10 kPa, the homogeneous broadening 7 « 
27r 860 MHz of the optical line exceeds the Doppler width kv ~ 27r 300 MHz, so velocity-changing collisions are 
inconsequential. The coUisional depolarization rate 7^ ~ 2n 70 MHz j20j is large compared to the inverse radiative 
lifetime I/te — 27r5.3MHz. The Rabi frequency \dE\/h k, I/te for the field intensity 8.8mW/cm^, which results in a 
saturation parameter {\dE\/K)'^ Tell ~ 10~^. The two-photon detuning is scanned in the range \Sb\ < 27r 1 MHz, and 
the ground-state relaxation rate can be estimated to be F « 27r 53 Hz ,3 1^ • 

Eliminating optical coherences with these approximations (for details see the Appendix) , we arrive at the following 
set of equations for the ground-state density submatrix {ogg = HgtrHg): 



(6) 



Tr{agJ = 1 , (7) 

where Ilg — X)m(|li (^i "^1 + |2, W (2, w|) is the ground-state projector, Ug = 2(2/ -I- 1) is the total number of 
sub-states in the ground state, and TTg is the total population of the excited state. The first term (oc tTe) of the source 
in © corresponds to the isotropic repopulation of the ground-state sublevels due to the spontaneous decay of the 
excited states. The other term (oc F) describes the entrance of unpolarized atoms due to diffusion and collisions. 
Due to the conservation of the total number of particles JT)), separate dynamic equations for the excited-state density 
matrix elements are not needed. Both the dynamics and steady state are completely governed by the non-Hermitian 
ground-state Hamiltonian: 

^cff = -^^(|l,m)(l,m|~|2,m)(2,m|) + i?-i^ng . (8) 



Here the excitation matrix, 



R= V |,^^/^,H(d-E.)^|e)(e|(d.E,)b-,^/)^^.^,^ 

^ h [{Sl - We) + 17/2] 



i,j ,e,m,m' 



E 



(»,m|(d.E,-)V)(e|(d.E,)|»,mO ^^^^^^ 
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contains the resonant (first summation) as well as off-resonant (second summation) contributions to the optical shifts 
and optical pumping rates (Hermitian and anti-Hermitian parts, respectively). The non-diagonal (i ^ j) elements of 
the resonant term induce the Raman coherence between the HF levels of the ground state responsible for the dark 
resonance. 

The generic matrix element in ^ is calculated from the Wigncr-Eckart theorem: 
{i,mMi^Et)'\e){e\{d^E,)lj.m,) = l{J,\\d\lJs,)l''-{F,,F,)r(F,.Fj) x 

X B-l)'^*'-"' { h F, } ^/SrFT(-l)'^- ( _ '< 5^ ) IE,- « E,)., , (10) 

where (Je||d||Jg) is the reduced matrix element of the dipole moment and 



r{Fe,F,) ^ V(2Je + l)(2^^e + l)(2F, + l) 

is the partial coupling amplitude of the Fi F^ transition. In the general case we have scalar {K = 0), vector 
{K = 1) and quadrupole {K = 2) contributions. All possible selection rules are contained in the coefficients of vector 
coupling, i.e., the 6j and 3jm symbols. 

For an atom moving along the direction of propagation of the optical field, the field frequencies are shifted due to 
the Doppler effect: uji ^ uji — kiV. As a result, a Doppler shift of the one-photon detuning 5l ^ ^l — kv occurs, where 
k = {ki + fc2)/2, as does a residual Doppler shift of the Raman detuning Sr — > 5p — (k^ — ki)v. At high buffer-gas 
pressure the residual Doppler shift is suppressed due to the Lamb-Dicke effect i2,^7\. However, in the general case 
the Doppler shift of the one-photon detuning can be significant, and certain quantities must be averaged over the 
Maxwell velocity distribution. Nevertheless, for buffer-gas pressures typically used in experiments, the approximation 
of homogeneous broadening is reasonable as a first approach to the problem, because the homogeneous width 7 equals 
or even exceeds the Doppler width kv. 

Here we consider the steady-state regime, setting {d/dt)dgg = in ((BJ. As a spectroscopic signal, we consider the 
total excited-state population TTe, which is proportional to the total light absorption in optically thin media or to 
the total fluorescence. The following procedure is used to find TTg. From lO, the ground-state density matrix dgg 
is expressed in terms of tTe, and then TTe is calculated from the normalization condition j?)). The solution of this 
algebraic problem can be obtained in a compact analytical form in two important special cases. The first arises when 
both field components have the same simple (circular or linear) polarization and there is no magnetic field. Here, for 
a suitable choice of the quantization axis, the excitation matrix R contains only diagonal elements with respect to 
the magnetic quantum number, i.e., m = m' in Q. The second case appears when a magnetic field is applied and 
just a few substates contribute to the Raman coherence for arbitrary light polarizations and arbitrary magnetic field 
directions. Both cases are considered below. 



\ Jg Je l\ 
\Fe F, I j 



III. SIMPLE LIGHT POLARIZATION, NO MAGNETIC FIELD 

We turn now to the case of circular field polarization when the quantization axis is directed orthogonal to the 
polarization vector (or alternatively linear polarization when the quantization axis is aligned along the polarization 
vector). We evaluate the total excited-state population, TTe, in order to determine how the dark resonance signal 
(proportional to TTg) depends on parameters such as the optical detuning from resonance. Under these assumptions, 
the complete set of equations (O can be split into independent blocks for each magnetic quantum number m (m- 
blocks). These blocks for m = ±Fi contain only one equation for the sub-state population tt^^^^\ The other blocks 
with m ±i^i contain four equations (two for the populations and two for the Raman coherences), corresponding to an 
effective two-level system with the upper |1, m) and lower |2, to) states (Fig.O. The parameters of the two-level system 
are expressed in terms of matrix elements of R as follows: the population relaxation rates = F + i?^^™'' include the 
optical pumping rates — 2 Im{(i, TO|i?|i, to)}; the dephasing rate is F12 — (Fi+F2)/2; the effective detuning 6ji — 
(5fl — (S'^'"'' — 5*2™'') includes optical shifts 5*1 = Re{(i, m|i?|i, m)}; and the coherence between levels is excited by the 
complex coupling t/ — i[/ = (1, m|i?|2, to). Note that the phase of the matrix element (l,TO|(d-Ei) |e)(e|(d-E2)|2, m) 
can be chosen equal to zero without loss of generality, so that (2,TO|i?|l,TO) = (1, TO|i?|2, m). 

Both the upper and lower states are repopulated with the same rate /3 — (TTe/re + F)/ng. First the total m- 

block population tt*^™^ ~ tTj™' + tTj™'' per unit repopulation rate is found. For the outermost blocks, to = ±fi, 
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FIG. 2: Effective two-level system, corresponding to one m-block 



detuning, 



where 



^i^^^O- The result for m ^ ±Fi is a quotient of polynomials of second order in the effective 



(m) _ 



B 



DSb + E 



(11) 



A = ri + r2 ; B = ri2[ri2(ri + r2)- 

c = rir2 ; Z? = 4C/y(fi - fa) ; 
E = r22fif2 + 2fi2(fi+f2)(F2_;72) 



8V^ 



(12) 



The repopulation rate, corresponding to unit total population in all m-blocks is 

-1 

-1*1 



.m=-Fi 



and the total excited-state population is finally expressed as 



Te{ngf3-T) . 



(13) 



(14) 



In the general case, when polarizations of the field components are different, or the same but elliptical, there is 
no basis where the matrices (1, m|i?|l, m'), (2, m|i?|2, m'), and {l,m\R\2,m') are simultaneously diagonal. In this 
situation, the full equation set for the ground-state density matrix elements must be solved, including all possible 
Zeeman and Raman coherences. Nevertheless, one important exception should be noted. If the optical linewidth 
is much greater than the excited-state HF splitting 7 ^ (wg^max — f^e.min), the quadrupole contributions to R are 
negligible 17]. The vector terms are diagonal (with respect to the magnetic quantum number) in the coordinate 
frame with z as the quantization axis, since [E^* x Ej] oc e^. Thus, we return to the case discussed above. 
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IV. DARK RESONANCES IN A MAGNETIC FIELD 

In a weak magnetic field, the ground-state magnetic sublevels are split due to the linear Zeeman effect, which can 
be described by the following additional term in the effective Hamiltonian (|SJ): 

Hb = ''^^'rnQ,i\i,m){i,m\ . (15) 

z,m 

Here the quantization axis is directed along the magnetic field, and Qi = fj,BgiB/h are the Zeeman splitting frequencies, 
with fj,B the Bohr magneton and B the magnetic flux density. The ^-factors of levels, 5^, are expressed through the 
electronic gj and nuclear gj Lande factors: 

.gj-gi , 

The magnetic field causes a precession of atomic coherences with frequencies m fl^ — m' ilj . When the Zeeman 
frequencies are much larger than off-diagonal elements of the excitation matrix f2i ^ m|i?|i, m')\, the light-induced 
Zeeman coherences within the i-th HF level are negligible. Thus, we again have a set of independent two-level systems, 
consisting of the sub-states \l,mi) and 1 2,7/12) (where |mi — 77/2] < 2 due to the selection rules). The formulas l|ll() 
and H12|l for the total block population are still valid for every {mi, m2)-block with the following substitutions: 

r. = r + Rt'^ ■ Sr = Sn - {s^r^ - st'^) {m,n, m^n,) ■ 

V - iC/ = (l,7ni|^|2,7772) = (2,7772|-R|l,7ni) . (16) 

If the Zeeman frequencies significantly exceed the widths F^, the Zeeman-split dark resonances are well resolved. 
In other words, the Raman coherence between the substates |l,7ni) and |2, 777.2) is effectively induced when the 
precession frequency is approximately equal to the Raman detuning: (5/? « nii Q.i — 70,2 ^^2- This condition can be 
simultaneously satisfied for only a few {mi, 77i2)-blocks. More precisely, the nuclear Lande factor is typically three 
orders of magnitude smaller than the electronic Lande factor (for cesium gj/gi ~ 2500); then, with good accuracy, 
ill = — ^^2 = ^ and the Zeeman shift of the dark resonance position is proportional to the sum of magnetic quantum 
numbers nfl = {mi -f 7/12) fi. It can be seen that, in the general case, three blocks (777, 777), (771 — 1, 7n -I- 1), and 
{m + 1, 777 — 1) contribute to the coherence preparation for the resonances with even shifts 2 m fl, and two other blocks 
(771 — 1, 77i) and {m, 7)7—1) contribute for the resonances with odd shifts (2 771 — 1) f2. When is tuned around the 
resonance with given shift 77 f2, the repopulation rate /3 can be written as 



-\ -1 



mi +7712— n 



where the first summand Z does not depend on the Raman detuning: 

^ 1 



, -n I -n I p(™2) / ' 

and 7r(™i'™2) jg (^j^g total population of the (7771, 7)72) block. Owing to the nuclear contribution, a further increase 
of the magnetic field causes the dark resonances to be eventually split into individual peaks, corresponding to each 
{mi, 7r72)-block [23l |. 

V. THE RESONANCE LINESHAPE 

We now consider the dark resonance lineshape in more detail. First, we analyze the particular case in which just 
two sub-states |1, 0) and |2, 0) participate in the Raman coherence, i.e., we consider the magnetically insensitive 
resonance (771 = 0) in a magnetic field. This (0, 0) resonance is of primary interest for possible clock applications 
because it is only sensitive to a magnetic field in second order. Here the absorption signal, 77^^^, 

form: 



1 F „ XT- / 1 1 



'■DR 



Z=yA——nr^+——n;r^ \ , (17) 



TeUg z + 7r(0) {5r) ng' ,f^o V r + ^l"' T + 4" 
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where tt^^-' is the total population of the (m = 0)-block per unit repopulation rate (see and (|12(l '). Since tt^"^ is 
a quotient of polynomials of second order in S^, the absorption can be written as the sum of an absorptive and a 
dispersive Lorentzian, and a constant background: 



(7/2)^ , ^ fe- '^0)7/2 

(7/2)2 + ~ 5nY ~^ ' (7/2)2 + (j^ _ 5^)2 



^Dfl = ,~/o^2 ^/ FT! + C2 ,~/o^2 , .7 \ ^2 + const . (18) 



The parameters in p8|) are expressed in terms of the coefficients introduced by in the following way. The dark 
resonance position is governed by the optical shifts and an additional term caused by the two-photon coupling between 
levels: 

The width of dark resonance reads 

The amplitudes of the symmetrical and antisymmetrical Lorentzians are found from the relations 

2 BC-AE-xAD 

^'^^'^^ = {A + czy ' (^^^ 

AD 

The background constant, C /{CZ + A) ~ T/ug, corresponds to the absorption far off the two-photon resonance. 

The result lO for the resonance lineshape is quite general. In fact, it does not depend on our simplified assumptions 
on the relaxation processes but is valid also in the low-saturation limit for arbitrary relaxation matrix, whenever only 
two states participate in the coherence preparation and Sji <C 7. 

Turning to the case of zero magnetic field and simple field polarization, we proceed with the goal of determining 
the resonance position, width and amplitudes of the symmetrical and asymmetrical components as above. Since all 
Zeeman levels within a given hyperfine level are now degenerate, we rewrite the repopulation rate /? H13(l as: 

(}^[z+{2F2 + l){n^'''\SR))^Y' ' (23) 

where 

Z = 




does not depend on 5r and corresponds to the absorption far off the two-photon resonance; the sum of the variable 
parts of the m-block populations tt^™' (5^) is expressed through the average over m-blocks, where the average of a 
variable X is defined as: 



Since n^"^\Sj^) is a quotient of polynomials of second order: 

(m) / r N "2 ^fl + "-^ 



n'-"'>{6R.) 



Sl + ai-^kR + bi"-^ 



im)_AD ^^m)_ BC-AE + ADis["''>~st^) 

"2 — -^yT ' "2 — Q2 ' y^^*) 

the average (tt'-^H'^-r))™ is a quotient of polynomials of order 2 {2F2 + 1). Generally this average describes a superpo- 
sition of resonances with different widths and positions due to the m-dependent power broadening and ac Stark shifts. 
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but if the laser detuning is not too large, \Sl\ < A, all resonances are well overlapped, and the average ( tt^™^ (J^j))^ 
can be approximated by a quotient of polynomials of second order. Here we use the following simple procedure, where 
the average of a quotient is substituted by a quotient of the averages: 

and where the correction factor a is chosen such that the exact and approximate expressions coincide at Sr = 0, i.e.. 



Our approximation for /3 yields an error less than a few percent across a wide range of parameters. With this 
approximation, we return to the resonance lineshape (|18|l . where the parameters are expressed in terms of the averages 
over m: 

(2F2 + 1)« (4™))™ 



So 

Z Zj A 

(7/2) Kb\ + (^2 >m - S^ 

Ci (7/2)2 ^ (2F2 + l)a>^2 )™ + (^2 )m'5o 



Z2 



C2(7/2) = (2F2 + l)a 



(a. 



2 



Z2 



1 r 

const = . (26) 



VI. COMPARISON WITH EXPERIMENT 



The analytical linesha pe (|18|l coincides exactly with the phenomenological model heuristically introduced previously 
to fit experimental data 18]. In those experiments a vertical-cavity surface-emitting laser (VCSEL) was modulated at 
the 9.2-GHz hyperfine splitting frequency of the cesium atom, so that the laser output spectrum contained modulation 
sidebands at this frequency. Using the carrier and one of the sidebands the dark resonance could be prepared and 
spectroscopically observed, as a function of the detuning Sj^ of the laser frequency from optical resonance. Data was 
taken for three different power ratios of carrier and sideband, with the cesium atoms contained in a cell with 8.7kPa 
of neon as a buffer gas. Detection used a modulation technique that allowed to extract simultaneously the absorption 
and the dispersion line shape For each detuning S^, both line shapes were simultaneously fitted by the model 
function H18|) . with Ci, C2, 7, and i5o as free parameters. Actually, as far as the line shapes themselves are concerned, 
this is a two-parameter fit: C2/C1 and 7 describe the shape, and the rest the overall amplitude and position of the 
dark line. 

Since these experimental data for Cs in Ne are fitted by (|18|l quite well, we can compare analytically calculated 
coefficients of the generalized Lorentzian to those extracted from experimental data. The dependence of the coefficients 
on the total light intensity X ol\E\^ ^\E2^ ''^^ almost trivial, at least when the power broadening [J^™^ -I- R^^'')/2 
exceeds the dephasing rate F in zero field: all the parameters Ci, C2, Sq, and 7 scale as I. Thus, the most representative 
test is provided by the dependence of the coefficients on the one-photon detuning 5^, and on the intensity ratio 
TZ = |i?i|2/|i52p between the two field components. Such comparisons with experimental fit parameters from 
are presented in Figs. El to El where Ci, C2, (5o, and 7 are plotted as functions oi 5l for three different relative 
intensities, TZ. The other parameters used in the calculations correspond to the experimental conditions: excitation 
by (7+ polarized radiation, total intensity I = 0.4mW/cm2, optical linewidth 7 = 27r750MHz, and ground-state 
relaxation rate F = 27rl50Hz. We use no free parameters, just a single trivial scaling factor for Ci and C2, and a 
constant offset for that accounts for the coUisional shift of the dark resonance position. 

We see a good qualitative agreement, especially for the resonance position Sq and for the width 7. There are some 
noticeable discrepancies for the amplitudes Ci and C2. In particular, we can see that the theoretical curve for Ci 
can cross the zero level at large Sl, which can be attributed to the well-known Raman absorption, but which is not 
observed in the experimental data. 
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FIG. 3: Absorptive coefficient, Ci, versus optical detuning, &l. Plots a), b) and c) are for TL — 2.4, 7.2, 22, respectively. The 
solid lines indicate the theoretical predictions while the points indicate the experimental data taken from 



VII. D2 LINE EXCITATION AND CONNECTION TO PREVIOUSLY EXISTING THEORIES 

In the specific case of the D2 line of Cs at high buffer-gas pressure, the two-photon amplitudes U and V are much 
smaller than the optical pumping rates i?^*^^ and the optical shifts 5*1"^ , respectively, because the most probable optical 
transitions Fi Fe = I+Je and F2 Fe = I — Je contribute to the one-photon transitions but not to the two-photon 
Raman coupling. Note that the ratio between V and R^^ can be arbitrary, depending on the one-photon detuning 
6l- As a result, the part of the absorption signal that varies with Sr is small compared to the constant one, and we 
arrive, to lowest orders, at the following approximate expressions. The parameter 

_ _D_ ^ 2(fi-f2)[/V 

2C- FiF^ 

is negligible with respect to the other contributions in (5o, 7 and Ci. The resonance position offset and the width are 
approximated as 

mr « ^«r?, + %^y^ (27) 
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FIG. 4: Dispersive coefficient, C2, versus optical detuning, 5l- Plots a), b) and c) are for IZ — 2.4, 7.2, 22, respectively. The 
solid lines indicate the theoretical predictions while the points indicate the experimental data. 



The amplitudes Ci and C2 are given by H21|l and 1)22(1 with a; = and 7 from H27|l. 

These results can be compared with those for a three-level A system in the low-saturation limit. Our formulas 
(ITHI)-(I22l will describe this last well, if we set Z = 0, i.e., 

C^mr^^^^; C,^/2=^. (28) 

Thus, the results are qualitatively similar (the main differences are the overestimated amplitudes Ci and C2), but 
now all parameters are unambiguously defined for the actual atomic structure. 

When C2 = the lineshape is symmetrical, and occurs if F = or Fi = F2. The first condition generalizes to 
Sl — 0, and the second corresponds to the condition of equal Rabi frequencies in a simple A system. 

When V — 0, the amplitude of the symmetrical signal is proportional to the square of the two-photon coupling: 

Ci « ^ J^V^^^ U' , (29) 



ri2 (Fi+Fa + ZFiF 



2 



•\2 



which is a key point of the perturbative studies jl7j but now, in addition, all effects of the optical pumping are 
accounted for in the prefactor in l|29|) . 
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FIG. 5: Frequency shift, Sq, versus optical detuning, 5l- Plots a), b) and c) are for TZ — 2.4, 7.2, 22, respectively. The solid 
lines indicate the theoretical predictions while the points indicate the experimental data. 



VIII. DARK RESONANCE POSITION. THREE POSSIBLE DEFINITIONS 



The center position of the dark resonance in essence determines the output frequency of the frequency reference or 
the magnetic field indicated by the magnetometer. Especially for asymmetrical resonances, it is somewhat unclear 
exactly how that center position is defined. The quantity Sq above is one possible definition of the resonance position, 
corresponding to the combined minimum of the absorptive part, and zero of the dispersive part, of the resonance 
described by (fTTji . 

Using H18|) - (|22|l . one can easily find another possible definition of the resonance center: the Raman detuning 
corresponding to minimum absorption. 

<5..n = - 4°^ + ^ . (30) 

A third possible definition is the point yo, where the dispersion associated with the absorption 118|l (by the 
Kramers-Kronig relations) is equal to zero. This is found to be 



(31) 
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FIG. 6: Dark resonance width, 7, versus optical detuning, Sl- Plots a), b) and c) are for TZ = 2.4, 7.2, 22, respectively. The 
solid lines indicate the theoretical predictions while the points indicate the experimental data. 



Each of these three quantities, Smin, Vo and Sq, could be considered the resonance center, depending on how the 
resonance is measured experimentally. In the general asymmetrical case, when V (non-zero effective one-photon 
detuning) and Fi ^ T2 (unbalanced optical pumping rates), all three values are different. Even their behavior versus 
Sl are qualitatively different (Fig. [Tj): near the one-photon resonance {V = 0) the centroid Sq of the Lorentzians 
has a dispersion-like shape, while Jmin is rather of an absorptive nature, and yo has a more complicated shape of 
mixed type. In addition, 60 and (5min are always finite, whereas yo goes to infinity at the zeros of Ci. These different 
dependences on optical detuning could, for example, alter the sensitivity of the frequency reference or magnetometer 
to the optical lock point. As a result, careful consideration must be given to the resonance detection method when 
designing frequency references or magnetometers based on dark resonances. 



IX. CONCLUSION 



Using very simple assumptions about the relaxation processes, analytical results can be obtained for the nonlinear 
absorption of bichromatic radiation near a two-photon resonance. The theory fully takes into account both the HF 
and the Zeeman level structures of alkali atoms as well as all light-induced effects. Our results constitute a good basis 
for understanding experimental works and further possible refinements of theory are possible. In particular, the case 
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FIG. 7: Three possible definitions of the dark resonance position. The centroid So corresponds to the sohd Une, (Smin - dotted, 
and yo - dash-dotted hne. All curves are calculated for the Cs D2 line. The parameters are T — 45piW/cm^, TZ = 0.5, and 
7 = 27r850MHz. 



of large Doppler width few > 7 can be immediately studied by the substitution 6^ 6^ — kv followed by averaging 
over the Maxwell distribution. 

In addition, the theory allows for a simple parameterization of experimentally measured dark resonances in terms 
of absorptive and dispersive components. The theory can therefore predict, for example, the detuning for which the 
dispersive part of the resonance is minimized and, for a given detuning, the asymmetry in the resonance lineshape 
that might be expected. The analysis of the different definitions of the resonance center position is also of interest for 
practical applications based on dark resonances such as atomic frequency standards and magnetometers. It appears 
likely that the additional understanding gained by the thorough theoretical analysis presented here will lead to further 
refinement and development of current and future applications based on dark resonances. 
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APPENDIX A: DERIVATION OF EQUATION © 

In this appendix we consider in detail the derivation of the basic equation set As is well-known the atomic 
density matrix obeys to the generalized optical Bloch equation. According to this equation, the evolution of the 
density matrix can be split into the two parts. The reversible one {d/dta = —i/h[H,a]) is governed by the total 
Hamiltonian of an atom in a resonant external field H = Hq + Hjj-e- The irreversible part originated from the 
interaction with environments (e.g. buffer gas or vacuum modes of electromagnetic field) are modeled by relaxation 
(super)operators of various kinds. The concrete form of the relaxation terms will be specified in the course of the 
derivation. ^ ^ 

The first stage is the elimination of the optical coherences a eg = IleCTlIg, where the operator He — 
l-Fg, TOg) (-Fe, 'Tiel projccts on the given HF component of the excited state. In the low-saturation limit the 
optical coherence matrix obeys the following equation in the rotating frame: 



-+7/2-z(<5i-We) 



'eg 



h 



ne(d • E,)n, 



^ne(d-E,)nje- 



'gg 



(Al) 
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On the left-hand side, the Raman detuning Sn is small compared to the homogeneous width 7 {\5r\ <C 7); Hi = 
X^m l^ii so that Ilg = III + 112. As is explained in the main text, the oscillations of the ground-state 

density submatrix agg can also be safely neglected in the rotating frame. Then, in the stationary regime (7t ^ 1) the 
solution of the equation ljAl|l is 



ne(d-E,)n» , V- ne(d-E,)nje-*("'-"^)* 



'eg 



^ I ■~r'n - ^('^L - t^e) ^ 7/2 - i((5L - We) - 




(A2) 



Under the conditions considered here, the equation for the ground-state density submatrix can be written 

jdgg = -r {dgg - ^^^^ ) - ^[Ho,dgg] " ^ (llg " h.c) + , (A3) 

where the line over operators indicates time averaging, i.e. all the oscillating terms should be removed from the 
product Hd~e^- Using ()A2f) . one finds that 



Ug Ho-E^^g ^ R^gg 



where R is the excitation matrix given by lO- The first term on the right-hand side of l|A3|l describes the relaxation 
in the ground state (due to both diffusion and collisions) toward the equilibrium distribution outside the laser beam, 

CTgg'' = Tlg/ug. All thc llucar (with respect to CTgg) terms, containing F, Hq, and i?, can be combined in the effective 
non-Hermitian Hamiltonian (|SJ| . The last term on the right-hand side of ljA3|l corresponds to the spontaneous radiative 
transfer of atoms from the excited-states, given by the density submatrix aee = Ve^Ve (where Ve — He), to the 
ground-state levels. Its structure will be specified below. 

In the low-saturation limit, the matrix aee obeys the equation 

d 1 i ^ ^ i / ^ \ 

■^^ee^-—^ee-^[He,^ee]-G{^ee}-^[VeHD-E^'Pe-h.C.j , (A4) 

where the first three terms on the right-hand side describe the radiative decay, the HF splitting {He — h ujeTle), 
and the coUisional depolarization of the excited state, respectively; the last term corresponds to the excitation due 
to light-induced transition from the ground-state levels. This last term can be considered as a source, because it is 
proportional to agg: 

1 (^ne'{d-E,)n,dggfij{d-Bj)^ne ^ 



h 7/2 + z(^L -CJe) 



E 



ne'(d-EOnjaggn,(d-E,)^ne 



7/2 + i((5L - t^e + tJi - Wj) 

The structure of the coUisional term G{^ee} can be found in Here we simply recall that during the course of a 
collision only the electronic component of the atomic polarization is depolarized. The nuclear component is involved 
in the process of depolarization due to the HF coupling. For all alkali atoms, the excited-state HF splitting Ag is 
much greater than radiative decay rate l/rg. In addition we assume that the coUisional relaxation rates 7^ for the 
excited-state electronic multipole moments of rank k = 1, . . . , 2Je -I- 1 also obey the conditions "f^Te S> 1 (for k = we 
assume 70 = 0, i.e. the collision-induced transitions between the fine structure components are not considered here). 
In this limit, AgTe ^ 1 and ■^K.Te ^ 1, the steady-state solution of HA4p has particularly simple form: 

f?ee = TTe — ; TTg = Te f 1 Tr{i? (Tgg } -|- C.C. j , (A5) 
Tie ^ ' 

which corresponds to total coUisional depolarization of the excited state. 

Here we shall illustrate this fact in one specific case, when the excited-state HF splitting is much larger than the 
depolarization rates 7k and when all the depolarization rates (except for 70) are the same (so-called pure electronic 
randomization model |23). If Ag » 7^, l/^e, one can neglect HF coherence in the excited state. For pure electronic 
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randomization both eigenvalues and eigenvectors of the LiouviUian Q are well-known |2(]| . which allows us to write 
the steady-state solution of 1A4|I for arbitrary JkT^: 



7k Te 



^ ^.l+lT.T^ ' (2Je + l) 



/ Fe Fe 


^ \ 




K 


\ 


I I I 


Je J 


{1 


I 


Je J 



nM{FeF,)Tr{fl^,{F^Fi)Se} . 



Here the source has the form 



(A6) 



ne(d-E0n,gg3nj(d-Ej) 

(7/2)2 + _ 



E 



ne(d-E0njgg3nj(d-E,) 

(7/2)2 _^ - + - ojjf 



the relaxation rates 



7L = 7k 



1-E 



(2Fe + 1)2 \ F, F, L 



(2Je + 1) 



I I Je 



i = 0, . . . , 2/ -I- 1 



(A7) 



correspond to the Zeeman projections of the nuclear multipole moments of rank L ,2Q] : and the Wigner tensorial 
operators are defined as 



TLAiiFaFb)^ J2 \Fa,ma)V2LTTi-r 



Fa -in a 



Fa L Fb 

—nia M mf, 



{Fb,mb\ 



As is seen from (|A7ll the rates 7^ are of the order of 7^ apart from 70 = 0. Then in the limit 7^X6 ^ 1 the leading 
term of (|A6(I corresponds to the summand with L ~ 0, which leads directly to the solution ljA5p . 

When the excited-state HF coherence is negligible, the radiative repopulation term in ljA3p can be written as 



A{^ee} = - ^ 



1 ^ r{F,,F,y 



TUF,Fi)deeTiq{FeFi) 



Fc,i,q 



One can easily prove the fundamental property: 



A{Ve} = --%. 

Te Ug 



(A8) 



(A9) 



which expresses the isotropy of the radiative relaxation. 

Thus, we see that in the case of total collisional depolarization of the excited state, when the excited-state density 
matrix is proportional to Ve [as shown in (|A5|) ]. ljA3p is reduced to lO- In addition, the expression for the optical 
coherence matrix (IA2|) allows one to calculate various spectroscopic signals (as well as the total absorption), for 
example, the total dispersion. 
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